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ABSTRACT 


We  present  numerical  results  obtained  with  the  mean 
trajectory  approximation  applied  to  a  non-adiabatic  process 
described  by  a  curve  crossing  model.  We  compare  them  to  the 
results  given  by  the  customary  trajectory  approximation  which 
uses  one  adiabatic  potential  to  g^enerate  the  nuclear  motion,  and 
to  the  results  given  by  Nikitin's  formula.  We  find  that  at  low 
kinetic  energy  the  three ’ procedures  give  very  different  results. 
The  differences  are  caused  by  the  fact  that  the  mean  trajectory 
method  includes  quantum  effects  in  the  effective  potential  used 
for  the  classical  motion.  We  also  discover  an  unexpected 
behavior:  under  certain  conditions  the  transition  probability 

seems  to  become  a  random  function  of  the  incident  energy. 


INTRODUCTION 


I  . 


The  description  of  the  dynamics  of  a  large  number  of 

surface  science  phenomena  requires  non-adiabatic  models  in  which 

the  system  undergoes  transltlonAbetween  at  least  two  states.  It 

is  customary  to  treat  such  problems  by  using  quantum  theory  for 

the  non-adiabatic  states  and  classical  equations  for  the  other 

T 

degrees  of  freedom.  If  the  energy  of  the  classical  variables  is 
low  compared  to  the  energy  difference  between  the  states,  this 
method  becomes  ambiguous  and  inadequate.  In  a  previous  paper  we 
reviewed  and  explained  this  situation  and  proposed  a  better,  but 
still  inexact,  method  which  we  call  the  mean  trajectory 
approximation.  In  this  article  we  present  a  deta.lled  numerical 
study  of  this  method. 

The  application  is  made  to  a  model  problem  in  which  a 
particle  approaches  a  surface  and  undergoes  a  transition  from  one 
electronic  state  to  another.  Since  the  particle  is  heavy  we 
intend  to  treat  its  nuclear  motion  by  a  classical  method. 

However,  since  in  the  Interaction  region  the  electronic  state  of 
the  particle  is  a  coherent  superposition  of  the  two  electronic 
states  mentioned  above,  we  do  not  have  a  unique  prescription  for 
writing  the  force  acting  on  the  particle  in  the  classical 
equation  of  motion.  By  using  systematic  methods  to  take  the 
classical  limit  for  the  nuclear  motion  MTA  arrives  at  a 
prescription  in  which  the  "classical”  potential  is  a  certain 
quantum  mean  of  the  electronic  matrix  elements.  This  is  a  time 
-dependent  potential  which  varies  as  the  collision  changes  the 
-p^jcupation  of  the  two  states.  If  the  state  i  is  highly  occupied 
throughout  the  collision  the  effective  potential  resembles  the 
potential  energy  curve  of  that  state.  If  the  occupation 

changes  from  1  to  2  the  effective  potential  switches  its 
character  from  resembling  to  resembling  course 

Intermediate  situations  in  which  the  potential  is  an  evenly 
weighted  combination  of  and  ^^e  also  possible. 
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The  present  paper  explores  numerically  the. properties  of 
the  effective  potential,  of  the  resulting  nuclear  trajectory  and 
the  differences  between  MTA  end  the  usual  trajectory 
approximation  (TA) .  One  of  the  more  striking  results  of  such 
calculations  occurs  when  the  coupling  between  the  diabatic  states 
is  small  and  the  kinetic  energy  is  low.  Under  these  conditions 
the  mean  potential  can* become  larger  .than  the  classical  total 
energy  of  the  particle  and  temporarily  trap  the  particle  at  the 
surface.  When  this  happens  the  asymptotic  value  of  the 
probabilty  of  emerging  on  a  specified  state  seems  to  become  a 
discontinuous,  random  function  of  the  kinetic  energy. 


II.  NUMERICAL  RESULTS  FOR  TWO  STATES  CROSSING 
II . 1  Summary  of  the  Equations 

2 

Since  the  theory  was  presented  in  a  previous  article  we 
give  here  only  the  equations  used  in  the  numerical  calculations, 
to  the  extent  needed  to  clarify  the  notation.  We  use  the  wave 
f^lnction 

i|>(R,X;t)  =  c^(t)  «j(x,R(t))  +  CgCt)  (2.1) 

where  x  represents  the  coordinates  of  the  electrons,  R(t)  are  the 
positions  of  the  nuclei  at  time  t,  i  »  1,2  are  electronic 

diabatlc  states  and  c^{t)  are  complex  probability  amplit.’des. 

The  ortho-normality  relation 

/  (x,R(t))  ♦p(x,R(t))  dx  «  (2.2) 

is  assumed  throughout.  The  amplitudes  c^(t)  satisfy  the 
equations 

ihc^  =  H^^(R(t))  c^(t)  +  H^p(R(t))  c^(t)  (2.3J 

with  ec=lor2,  8  =  1  or  2  and  p  p:  a,  and  c  =  dc  /dt. 

a  a 

For  the  matrix  elements  H^p{R(t))  we  chose 

H^^{R(t))  -  U^{(Rg/R)  expC-2{R-R^)]  -2exp[-(R-R^) ] }  +  X/2 

(2.4a) 

H22(R(t))  »  U^{(R^/R)  exp[-2(R-R^)]  +  2expt-(R-R^) ] )  -  x/2 


(2.4b) 


(2.4c) 


H,,(R(t))  =  H-,{R(t))  =  -2pU  exp[-(R-R^)] 

12  2  1  O  0 

The  functions  Hjj(R)  and  are  the  diabatic  energy  curves  of 

the  two  states  #^(x.R(t)).  They  are  plotted  in  Fig.  1  as  full 
lines.  The  curve  is  a  Morse  potential  (which  we  aultlplied 

by  (R  /R).  to  nake  the  repulsive  energy  infinite  at  short  surface 

O 

particle  distances  R)  to  which  we  added  a  constant  energy  Ae/2. 
The  curve  is  modified  to  have  smaller  binding  energy  than 

Hjj.  The  coupling  dies  off  exponentially  on  the  same  len’gth 

scale  as  the  potentials  and  becomes  constant  at  short  distances. 
The  coupling  constant  is  p. 

The  equations  (2.3)  and  (2.4)  are  written  in  terms  of 

dimensionless  variables.  The  length  unit  a~^  is  given  by  the 

inverse  length  os  appearing  in  the  Morse  potential  V(y)  >• 

U^{exp[-2a(y-y^)]  -  exp ( -ot( y-y^ ) ) ) .  The  length  o”^  is  the  range 
o  o  o 

of  the  potentials  H,,  and  H*..  The  time  unit  is  (osv  )”  ,  where 

22  ■  0 

v^  is  the  velocity  of  the  incident  particle.  For  energy  we  use 

hosv  and  for  mass  (ha/v  ).  An  important  dimensionless  quantity 
o  0 

is  the  Massey  parameter 

X  -  Ae/hav^  ,  (2.5)- 

where  Ae  is  the  asymptotic  energy  gap  between  the  states  (i.e.  Ae 
=  lim(Hj^j^(R)  -  H22(R))  for  R  **■  *)  •  The  values  of  the  parameters 
used  in  the  present  calculations  are  given  in  Table  I. 

We  find  it  useful  sometimes  to  switch  between  the  diabatic 
"and  the  adiabatic  pictures.  Therefore  we  present  briefly  the 
Equations  used  to  make  this  connection.  Since  the  time 
dependence  appears  through  R  only  we  can  write  the  Eqs .  (2.3) 

(using  an  obvious,  two  dimensional  matrix  notation)  as 

(  ihv)  3C/3R  <•  T*  ‘  "c  (2.6) 

with  V  H  R.  Since  we  have  assumed  that  the  diabatic  states  are 


orthonorm&l  we  can  use  an  unitary  matrix 


U 


wi  th 


'  cose  ;  sine  ' 
.-sine  ;  cose 


tan  2e 


2H, , 
12  22  11 


The  adiabatic  surfaces  are  given  by 


(2.7) 


(2.8) 


(R)  -(U“^)  H  U  ,  a  -  1  or  2,  (2.9) 

eta  '  ay  yH  |io:  \  ^  i 

(repeated  Indeces  are  summed  over)  and  the  adiabatic  amplitudes 
by 


‘^a  a  -  1  or  2  .  (2,10) 

Defining 

R 

c“  5  C^  exp{(i/h)  J-  (H®^{R)/R)dR)  (2.11) 

0* 

and  using  the  above  equations  in  (Eq.2.3)  leads  to 

R 

3c  /3R  =  (-l)“(3e/3R)  c„  exp{(l/h)  J  (H®  -H®  ]  R'^dR)  . 
o  ^  ^  p^  aa 

(2.12) 

where  a  1  or  2,  p  1  or  2  and  p  a. 

In  order  to  solve  these  equations  in  either  representation 
we  must  specify  R(t).  The  mean  trajectory  approximation  (MTA) 
prescribes  the  use  of  Hamilton  (or  Newton)  equations  with  the 
effective  potential 


V(  R  ) 


(2.13) 


i  .  j  =  .1 


In  what  follows  we  solve  siaultaneously  the  equation  of 

for  R(t)  and  that  for  the  aaplltudes  c  (or  c  ,  In  the 

ce  a 

repr es en ta t i op  )  . 


motion 

adiabatic 


I I. 2.  NUMERICAL  RESULTS 


We  present  here  the  mean  trajectories,  the  occupation 
probabilities  of  the  states  1  and  2.  and  the  effective  potently’ 
as  functions  of  time,  for  various  values  of  the  Massey  parameter 
X.  The  results  are  compared  to  Nikitin's  formula^  and  to  those 
obtained  by  running  trajectories  on  the  lowest  adiabatic  state. 

The  relationship  between  the  X  values  used  in  these 
calculations  and  the  incident  kinetic  energy  is  given  in  Table 
2.  The  mass  of  the  Na  atom  and  the  parameter  values  indicated 
in  Table  1  were  used  to  make  the  connection.  We  have  mv  /2  - 
(nAe^/2ti^a^X^)  =  3.30  X“^  (Hartree)  -  89.9  x"^  eV. 

The  numerical  calculations  consist  of  solving  the 
differential  -equations  (2.3)  (for  the  diabatic  representation)  or 
(2.12)  (for  the  adiabatic  representation),  simultaneously  with 
Newton's  equation  with  the  effective  potential  Eq.  (2.13). 

The  mean  trajectory  R(t)  for  a  particle  starting  in  state  1 
is  shown  in  Fig.  2.  From  each  curve  we  can  obtain  the  values  of 
the  following  parameters:  the  position  and  the  tine  at  which  the 
particle  turns  (R  changes  sign)  and  the  times  when  the  particle 
goes  for  the  first  and  the  second  time  through  the  diabatic 
crossing  point  R  =  5.  All  these  parameters  characterizing  the 
kinematics  on  the  mean  trajectory  show  Interesting  quantum 
behavior.  If  we  were  to  run  a  Newtonian  trajectory  on  the  lower 
adiabatic  state  we  would  expect  all  the  trajectories  considered 
‘.here  to  turn  at  approximately  the  same  distance  R.  located  on  the 
sharp  repulsive  part  of  the  surface.  Fig.  1  places  this  value  of 
R  between  1.2  and  1.6.  In  Fig.  2  we  see  that  the  value  of  R  at 
the  turning^  point  has  the  values  3.3.  2.15.  1.8  and  0.7  (for  X  » 
20.  10.  5  and  2,  respectively).  Furthermore,  the  times  of  the 
second  crossing  are  7.6.  10.  11.2  and  12.3  for  X  -  20,  10,  5,  2. 
respectively);  they  have  the  striking  feature  that  the  low  energy 


particle  (e.g.  X  =  20)  recrosses  much  earlier  than  the  fast  one 
(e.g.  X  -  2) . 

•  The  two  effects  are  related,  since  the  low  energy  particle 
reaches  the  second  crossing  first  partly  because  it  turns  around 
earlier  than  the  high  energy  one.  These  observations  can  both  be 
rationalized  as  quantua  effects  introduced  through  the  effective 
potehtial  which  depends,  through  the  electronic  aaplltudes,  on 
the  incident  energy  (hence,  on  X).  As  shown  in  Figs.  (3.a)- 
(3.d)  the  effective  potential  can  differ  markedly  from  the 
adiabatic  ones.  Since  we  start  with  a  particle  in  the  diabatic 
state  1  the  effective  potential  at  large  R  resembles  .  As  the 
particle  approaches  the  surface  and  the  transition  1  -►  .2  starts 
taking  place,  the  effective  potential  differs  from  both  adiabatic 
ones.  If  the  population  of  the  adiabatic  state  2  is  larger  than 
that  of  the  state  1  the  effective  potential  resembles  ..  This 
is  roughly  the  case  in  Fig  3.  where  X  -  2  corresponds  to  a 
kinetic  energy  of  22.4  eV.  Lowering  the  kinetic  energy  gets  the 
effective  potential  closer  to  H® ^  .  Since  the  effective  potential 
can  be  computed  only  at  points  reached  by  trajectories  the  curve 
stops  at  the  values  of  R  corresponding  to  the  turning  point 
(given  in  Fig.  2).  Finally,  note  that  the  effective  potential 
for  the  incoming  trajectory  is  different  from  the  outgoing  one; 
this  happens  because  the  amplitudes  c^  are  different  on  the  two 
trajectories . 

The  reason  why  both  the  turning  and  the  recrossing  times  go 

-down  with  energy  can  be  understood  by  looking  at  Figs  (3. a)  - 

"(3. d).  At  the  highest  energy  V  is  close  to  H®  ,  while  at  low 

^ *  2  2 

"energy  V  is  close  to  H, ,  .  Thus,  an  energetic  particle  has  to 
e  r  I  ji  1 

travel  a  longer  distance,  both  to  turn  around  and  recross,  than  a 
slow  one . 

The  probability  of  remaining  on  the  Initial  adiabatic 


state  is  plotted  in  Fig.  '4.  for  four  values  of  X.  All  the 
curves  change  rather  suddenly  as  the  particle  passes  through  the 
crossing  region  (R  -  5).  The  magnitude  of  this  change  increases 
monotonical ly  with  X  that  is  with  the  kinetic  energy.  The 
largest  kinetic  energy  gives  the  biggest  .drop  in  P*.  This 
correlates  nicely  with  the  fact  that  in  Pigs.  {3. a)  -  ( 3 .  d)  the 

effective  potential  which  starts  by  being  close  to  , 

*  21  •  • 

crosses  over  to  resemble  snail  values  of  X,  but  renains 

similar  to  for  X  ■  10  .or  20.  The  behavior  of  P®  at  the 

second  crossing  is  not  a  monotonic  function  of  X.  For  two  values 

of  X  (i.e.  20  and  5)  P*  decreases  at  the  second  crossing,  while 

for  X  equal  to  10  or  2.  P^  increases  at  the  second  crossing. 

This  suggests  that  the  asymptotic  value  (i.e.  the  post  collision 

value)  of  P^  oscillates  with  X.  That  this  is  the  case  is  seen  in 

Fig.  5  where  we  plot  the  final  values  of  P®  as  a  function  of  X, 

for  values  of  X  around  X  •  2.1  (i.e.  incident  kinetic  energy 

around  20.3  eV).  On  the  same  graph  we  plot  the  results  of 

3 

Nikitin's  theory  which  uses 


P  -  2P,^2(1-P^^2). 


exp(-X(l-r))  (sinh[X( l+r) ]/  sinh[2X]> 


(2.14) 

(2.15) 


and 


,,  „2,-l/2 

7  -  (  1  ♦  p  ) 


(2.16) 


The  Nikitin  formula  does  not  give  rise  to  oscillations,  but  it 
goes  smoothly  through  them.  Reasonable  qualitative  agreement 
with  our  numerical  results  is  obtained  if  Nikitin's- equation  f.or 
P  i-s  multipi  ifrd  ’by  V.  ’  ‘ 


t. 


2sin‘[  ; 

t , 


(R(t))  -  (R(t))]dt/2h  +  4>]  (2.17) 


where  0  is  an  adjustable,  time  independent  phase,  t,  is  a  pre 


3 


collision  time  and  is  the  timfr  at  which  the  trajectory  reaches 

the  turning  point;  this  factor  is  suggested  by  arguments  similar 

to  chose  made  in  the  classic  Stuckelberg  paper  and  are  not 
4  ^ 

repeated  here. 


In  order  to  find  the  period  of  oscillation  we  use  dt  » 
dR/lv|,  where  v  »  R,  and  the  expressions  for  and  11^2  given  by 

the  Eqs.  (2.4).  Assuaing  that  |v|  is  constant  we  obtain 


R 

I  s  : 

R, 


(Mvl)'^  s  {4U  exp(-a(R-R  )  -  Ae)^  +  exp(-2a(R-R  )}^^^dR 

o  o  O  o  .0 

(2.18) 

Here  I  is  the  tine  dependent  tern  in  the  argument  of  the  sin 
function  in  the  previous  equation.  Since  the  major  contribution 
to  this  integral  cones  from  values  of  R  near  the  turning  point, 
we  neglect  Ae .  This  allows  us  to  perform  the  integral.  Since  R^ 
is  arbitrarily  large  we  chose  it  so  that  all  the  terms  exp  (-o:Rj^) 
are  zero.  After  all  these  we  obtain 

I  -  -20  (1*8^)^'^^  (X/Ae)  exp(aR  ) 
o  o 

2 

Thus  sin  {I*<^)  oscillates  as  a  function  of  X  with  a  period  of 
2vX/I  (this  is  independent  of  X!).  The  parameter  values  used  in 
the  present  work  lead,  when  Inserted  in  the  above  equation  to  a 
period  of  0.081  which  is  close  to  that  given  by  the  numerical 
calculations.  This  might  seen  gratifying,  since  such 
oscillations  are  in  principle  observable  and  the  simple  analysis 
presented  above  links  their  period  to  the  asymptotic  gap  Ae,  the 
binding  energy  of  the  ion  ,  the  strength  of  the  coupling  p.  and 
the  range  of  the  potential  a  Unfortunately  in  most  surface 


science  problems  the  phase  will  be  randomized  to  soms  extent  by 
inelastic  processes  and  this  will  smooth  out,  of  perhaps  even 
wipe  out,  the  oscillations. 

The  behavior  of  a  system  which  starts  in  the 'adiabatic 
state  2  is  shown  in  Figures  6,  7. a  -  7. d  and  8.  The  discussion 
of  these  results  follows  along  the  lines  of  the  case  when  the 
particle  starts  in  the  state  1,  and  it  is  therefore  omitted. 

A  comparison  between  the  mean  trajectory  approximation  and 
a  theory  in  which  the  particle  is  propagated  on  the  adiabatic 
surface  2  (the  Initial  state  is  2)  is  presented  in  Figs.  9. a  - 
9.d.  The  fact  that  the  largest  differences  between  'the  two 
theories  are  obtained  for  the  smaller  values  of  X  (i.e.  X  «  2  and 
5)  might  seem  surprising,  since,  in  general,  one  expects  the  two 
theories  to  differ  less  at  high  energy.  However,  an  examination 
of  the  Figures  7.  a  -  7.d  shows  the  reason  for  this*  result:  the 
effective  potential  at  small  values  of  X  resembles  fairly 
closely,  and  differs  from  ^^2'  small  X  the  mean 

trajectory  differs  most  from  a  classical  trajectory  on  ^^2' 
discrepancy  between  occurs  because  (Fig.  8)  goes 

down  at  the  first  crossing  for  X  ®  2  euid  X  =  5  much  more  than  for 
X  s  10  or  20.  Th\is  the  effective  surface  switches,  from  being 
close  to  Hgg  prior  to  the  first  crossing,  to  becoming  close  to 


If  we  attempt  to  apply  Nikitin's  formula  we  must  take  into 
account  the  fact  that  the  velocity  in  the  region  where  the  two 
states  interact  strongly  differs  from  the  incident  velocity.  We 
include  this  effect  by  using  Nikitin's  formula  with  X*  »  de/hav  , 

W 

instead  of  X  •  de/hav  ,  where  v  is  the  velocity  at  the  crossing 

®  2^2 

point.  This- is  given  by  mv  »  mv^  -  de  and  as  a  result  X*  can  be 
written  as  X*  »  X(l-6X^)  ,  with  6  =  (?ia)^,/mde.  The  transition 

probability  after  the  first  crossing  is  now  given  by 


(2.19) 


exp[-X  (l-r)l[sinh  X  {l+y)]/[slnh  2X  ] 


with  y  defined  by  (2.16).  In  Fig.  10  we  show  the  dependence  of 
^2-*i  ^  given  by  the  mean  trajectory  approximation  (broken 

lizxe)  and  by  Eg.  (2.19)  (full  line),  for  several  coupling 
strengths  ^.  The  quantities  and  X^  marked  on  the  graphs  are 
defined  by  X^  =  Ae/hcev^^  and  X^  ■  Ae/hocv^  where  v^^  corresponds  to 

a  kinetic  energy  equal  to  the  asymptotic  mismatch  Ae  (i.e.  v  « 

1/2  ^ 
(2Ae/m)  )  and  v^  to  a  kinetic  energy  equal  to  the  difference 

between  the  potential  energy  at  the  crossing  point  and  the 

aisymptotic  energy  of  the  lower  energy  curve  (i.e.  v-  =« 

1/2  ^ 

(Ae/m)  ).  Nikitin's  formula  for  ^2-*i  deviates  only  slightly 

from  the  mean  trajectory  values  as  X  exceeds  X^^ .  From  Nikitin 

formula  =  0  if  2:  5^^  while  the  mean  trajectory  calculations 

give  small,  finite  values  for  P  _  if  X  >  X_  which  go  to  zero  as 

X  is  increased. 


Nikitin's  formula  is  ftirther  improved  if  we  use  in  Eq.  2.19 

««  « 

X  =  Ae/hctv^  Instead  of  X  ,  where  v^  is  obtained  from 
P  P 

mVp/2  +  (Rp)  »  »v^/2  -  Ae/2  (2.20) 

and  R  is  the  position  where  3e(R)/3R  is  maximum.  For  the  model 
P  "  91/9 

xised  here  (i.e.  Eqs.  (2.4)  R^  =»  R^  +  a  ln[4(l+p  )  '  W^/Ae], 

with  R^  is  defined  by  (2.4).  These  lead  to 
•  •  *•2  —1/2 

X  »  X(1  -  6X^)  (2.21) 

'.and 

6  -  6[1  -  [2(l-y)]^''^  +  (r^/16)(Ae/0^)  (R^/Rp)]  .  (2.22) 

«•  • 

Using  Eq.  (2.19)  with  X  instead  of  X  gives  almost  perfect 
agreement  with  the  meem  trajectory  approximation. 


2-**l 


Unfortunately  P 


is  not  measurable  and  the  use  of 


Nikitin's  formula  to  compute  P®,  the  probability  that  the  system 
ends  up  in  the  adiabatic  state  1  after  having  started  in  state  1, 
can  be  done  with  moderate  success  only  after  a  correction  for 
interference  effects  are  made.  To  do  this  we  multiply  Nikitin's 
formula 


(2.23) 


with  the  factor  given  by  Eq.  (2.17)  which  includes  interference 
effects.  The  computation  of  uses  Eq.  (2.19)  with  X** 

instead  of  X  .  We  use  this  corrected  formula  for  a  qualitative 
discussion  of  our  numerical  results  and  for  madcing  crude 
estimates. 


In  Figs.  11  and  12  we  show  the  mean  trajectory,  results  for 
P^  (full  lines),  together  with  the  results  obtained  by  running  a 
trajectory  on  (broken  lines)  and  those  given  by  Eq.‘(2.23) 
uncorrected  for  interference  (i.e.  Eq.  (2.23)  without  multiplying 
it  with  (2.17).  We  see  that  as  the  coupling  constauit  p  and  the 
incident  kinetic  energies  are  decreased  the  results  obtained  with 
the  mean  trajectory  approximation  differ  markedly  from  those 
obtained  by  running  a  trajectory  on  •  This  happens  because, 
the  mean  potential  differs  substantially  from  .  Qualitatively 
this  can  be  understood  as  follows:  since  p  is  small  the  diabatic 
states,  are  weakly  coupled  and  the  system  tends  to  remain  on  its 


initial  diabatic  state.  In  the  present  case  the  system  starts  on 
^22  large  R  coincides  with  Hg,)  aiid  tends  to  stay  on 

it.  Since  for  small  values  of  R  ^  this  means  that  the 

system  tends  to  undergo  a  transition  (in  the  adiabatic  picture) 
from  to  ^11*  mean  potential,  which  is  initially  equal  to 

begins  to  resemble  more  than  H22*  Thus  after  the 
crossing  point  *»  and  therefore  the  calculation  using  H22 

to  obtain  the  trajectory  leads  to  results  which  are  very 
different  from  the  mean  trajectory  ones. 


In  the  remainder  of  this  section  we  investigate  a 
remarkable  phenomenon  that  takes  place  for  sball  p  values.  It  is 
clear  that  the  effective  potential  varies  in  time  as  driven  by 
the .  qucintvim  au^lit^des  c^.  Imagine  a  trajectory  whose  initial 
total  energy  is  close  to  the  maximum  of  *  transition  to 

the  adiabatic  state  1  occurs  the  effective  potential  changes  into 
a  mix  of  and  •  Since  >  E^^  this  can  have  the  effect 

of  making  the  barrier  of  larger  them  the  barrier  of  H22  smd 

the  total  initial  energy  of  the  particle.  It  is  thus  possible 
that  a  trajectory  that  would  pass  over  that  surface  was 

used  in  Newton’s  equation  ~  can  be  reflected  by  since 

that  potential  will  have  a  larger  barrier.  Furthermore,  a 
trajectory  which  passed  over  the  barrier  of  in  its  way  to 

the  surface  may  be  trapped  in  its  way  out  by  a  surge  in 
caused  by  an  Increase  of  the  transition  amplitude  to  the  upper 
adiabatic  state  1.  Due  to  the  all  or  nothing  character  of  the 
classical  reflection  process  we  expect  to  see  some  sudden  changes 
in  the  behavior  of  the  system  as  we  vary  p  (or  X)  through  the 
range  in  which  such  a  process  can  take  place.  The  numerical 
results  show  the  onset  of  a  remarkable  behavior:  within  our 
ability  to  resolve  the  issue  the  transition  probability  P® 
becomes  a  random,  discontinuous  function  of  X.  We  expect  the 
phenomena  mentioned  above  to  take  place  when  the  effective 
potential  exceeds  to  total  energy  of  the  particle.  After  the 
second  crossing  the  effective  potential  is  roughly 

1/2  AePf  -  (Ae/2)(1-P®)  =  Ae(P®-l/2)  and  the  asymptotic  energy 

^  ^  2  ■  2  a 

conservation  condition  is  bv^/2  ^  1/2  Ae  =  mVj/2  +  Ae(P^  -  1/2). 

It  is  possible  to  obtain  zero  final  velocity  Ae(P*  -  1/2)  mv^/2 

-  Ae/2.  This  leads  to 

pf  k  mv^/2Ae  =  (26X^)"^  .  :  - .  -  ( L  24 )"'  '  “ ' 

1  o 

This  equation  permits  us  to  estimate  under  what  conditions  we  can 
expect  the  effective  potential  to  turn  the  particle  back  towards 
the  surface  after  the  second  crossing. 


In  Fig.  13  we  plot  the  maximum  value  obtained  max  P®  of  P® 
obtained  from 

-  ^2-1^  '  (2.25) 

with  ^2-*l  Siven  by  Eq,  (2.19).  The  factor  of  4  appears  because 

the  Nikitin  formula  corrected  to  Include  interference  is  pf  *•  4 

2  ^ 

^2-*i  “  ^21^  ®  and  the  maximum  value  is  obtained  taking 

the  sin  equal  to  one.  The  broken  curve  is  (2£X^)’’^  and  the  full 

lines  are  the  maximum  values  of  computed  as  described  above, 

for  various  values  of  ^.  For  ^  0.1  the  broken 'and  the  full 

lines  cross  and  the  condition  (2.24)  is  satisfied  in  the  range  of 

X  values  indicated  by  arrows  (for  p  =  0.1).  This  is  roughly  the 

range  (X^^,  X^),  where  X^^  and  X^  were  defined  earlier.  For  higher 

values  of  coupling  constant  p  (2.24)  is  not  satisfied. 

While  we  have  used  approximate  formulae  to  obtain  the  plots 
in  Fig.  13  they  give  us  a  fairly  reliable  hint  concerning  the 
parameter  region  where  we  should  look  for  unusual  behavior  in  the 
mean  trajectory  results.  In  Fig.  14a  we  show  the  dependence  of 
P®.  given  by  the  mean  trajectory  approximation,  for  p  =  0.1  and  X 
*  30/  30.015  and  30.020.  P2(t)  is  the  probability  that  the 

particle  is  on  the  adiabatic  state  2  at  time  t.  The  horizontal 
broken  line  indicates  the  values  of  P®  below  which  the  condition 
(2.24)  is  satisfied.  For  X  30  and  30.015  the  condition  is  not 
satisfied  at  any  time  and  the  behavior  of  P®  is  normal:  at  t  --  7 
the  particle  undergoes  its  first  crossing  and  at  t  11  the 
second  one. 

Since'  the  asymptotic  value  of  P®  oscillates  as  a  function 
of  X  the  two  curves  go  to  different  asymptotic  values.  The  curve 
X  »  30.02  behaves  quite  differently.  After  the  second  crossing 
at  (t  11)  P®  satisfied  Eq.  (2.24)  and  the  effective  potential 


exceeds  the  total  energy  of  particle.  Therefore  the  particle, 
which  is  aoving  away  from  the  s\irface,  is  turned  back  and 
undergoes  a  third  crossing  at  t  -  38  which  takes  above  the 
broken  line;  the  fourth  crossing  lifts  even  higher  so  the 
energy  of  the  particle  exceeds  that  of  the  effective  potential 
and  the  particle  escapes.  This  interpretation  is  substaoitiated 
by  the  trajectory  plots  shown  in  Fig.  14b,  for  the  same 
conditions.  The  crossing  point  is  at  R  »  5  and  the  trajectory 
corresponding  to  X  =  30.02  crosses  it  four  times.  In  Fig.  14c  we 
show  the  effective  potential  for  X  «  30.02,  plotted  together  with 
the  adiabatic  potentials  (broken  lines) .  As  the  particle 
approaches  the  surface  the  effective  potential  is  the  curve 
indicated  by  (1).  After  the  crossing  point  R*5.0  the  potential 
has  the  form  shown  as  (2).  This  curve  stops  at  the  turning  point 
and  the  height  of  that  point  equals  the  total  energy  of  the 
particle  (since  there  the  velocity  is  zero).  The  potential  (2) 
is  half  way  between  and  it  is  repulsive.  After  the 

turning  point  the  probability  P®  does  not  change  until  the 
particle  reaches  crossing  point  and  the  potential  remains 
unchanged  (i.e.  it  is  given  by  (2)).  At  the  crossing  point  the 
probability  changes  and  the  effective  potential  (3)  resembles 
more  H®^  and  it  becomes  attractive.  At  R  a  bit  larger  than  7  the 
particle  reaches  a  turning  point  and  goes  back  towards  the  wall. 
At  the  crossing  the  probability  changes  again  and  the  effective 
potential  is  given  by  curve  ( 4 ) .  This  resembles  more  the 
potential  H®^  than  the  effective  potential  (2)  which  acted  on  the 
particle  when  it  first  approached  the  wall.  The  particle  hits 
the  wall  and  turns  around  and  the  effective  potential  is  still 
.(4)  since  P®  is  not  changing.  At  the  turning  point  P®  becomes 
very  small  and  the  effective  potential  (5)  is  repulsive  (since  it 
is  largely  composed  of  Hjj)  and  the  particle  escapes. 

To  investigate  this  phenomenon  in  more  detail  by 
calculating  the  asymptotic  values  of  P®  as  a  function  of  X,  for  X 
in  the  range  [30.00.  30.2]  in  which  the  unusual  behavior  takes 


place.  The  results  ajrs  represented  In  Pig.  15.  In  Fig.  .15a  we 
shew  by  full  lines  the  mean  trajectory  results  for  the  case  when 
the  behavior  is  normal.  in  that  range  is  a  continuous 

function  of  X.  The  circles  represent  results  obtained  with  the 
mean  trajectory  approximation  in  the  range  of  X  values  where 
multiple  trajectories  occur.  The  numbers  near  the  circles 
indicate  how  many  times  the  trajectory  was  turned  back  towards 
the  wall.  If  the  number  is  1  the  trajectory  hit  the  wall  twice, 
if  the  number  is  4  then  the  the  wall  was  hit  five  times,  etc. 

The  broken  lines  show  the  results  for  P®  at  the  time  when  the 
particle  is  moving  away  from  the  surface  hits  the  wall  for  the 
first  time  and  hits  the  turning  point  that  will  send  it  back 
towards  the  surface  for  a  second  collision.  The  horizontal  line 
is  given  by  the  Nikitin  formula,  uncorrected  for  interference. 

In  order  to  find  whether  the  circles  follow  a  pattern  we  carried 
out  mean  trajectory  calculations  in  a  very  narrow  range  Xe[30.01 
30.05].  Even  with  this  very  high  resolution  we  still  find  that 
the  values  of  P^  represented  by  the  circles  seem  to  originate 
from  a  discontinuous  function  f(X)(l.e.  P®  =  f(X)).  In  Fig.  16 
we  show  that  the  phenomenon  appears  for  other  values  of  X. 


IV 


DISCtlSSibN 


The  mean  trajectory  approximation  introduces,  through  the 
mean  potential,  some  quantiim  effects  into  the  "classical" 
trajectory.  To  understand  how  such  effects  work  and  how  we  can 
improve  their  treatment  it  is  useful  to  recall  one  of  the 
derivations  of  the  MTA.  We  have  shown  that  we  can  generate  MTA 
by  running  a  narrow  wave  packet  on  the  two  state  systems.  This 
is  done  by  chosing  a  wave  function  of  the  form 

«|i(x,R;t)  =  G(R;  t)  {c^<t^(x,R(t)  )  +  c^'tj  {x,R(  t )  )  } 

where  G(R;t)  is  a  Gaussian  wave  packet.  Using  a  method  developed 
by  Heller^  we  can  generate  equations  for  the  motion  of  the  center 
of  the  wave  packet  which  turn  out  to  be  the  mean  trajectory 
eqfuation.  These  equations  will  attempt  to  mimic  the  quantum 
motion  of  the  packet.  The  temporary  trapping  of  the  particle  can 
be  rationalized  in  terms  of  this  picture  as  follows.  If  the 
packet  passes  the  crossing  point  and  P®  is  very  small,  it  will 
move  like  a  classical  particle  on  the  lower  adiabatic  state  ^^2' 
Let  us  assvime  that  the  particle  hits  the  surface  and  when  it 
cones  back  becomes  large.  The  packet  now  moves  mostly  on  the 

upper  state  and  it  is  trapped  in  its  potential  well.  As  it 

rattles  around  it  passes  through  the  crossing  point  and  if  Pj 
stays  high  the  packet  remains  trapped;  if  P^  becomes  small  the 
packet  goes  on  ^^2  escapes  (if  it  happens  to  move  away 

from  the  wall)  or  it  collides  with  the  wall  (if  it  happens  to 
move  towards  it) .  If  the  collision  with  the  wall  takes  place  the 
particle  turns  around,  it  may  repeat  the  above  sequence  (if  P^ 
becomes  large),  or  it  may  get  away  (if  P®  stays  small). 

Attempting  to  describe  the  motion  of  the  particle  by  using 
one  packet  agrees  in  a  superficial  way,  with  our  intuition.  Me 
want  to  treat  the  atom  notion  classically  and  it  is  well  known 
that  we  can  achieve  that  by  taking’s  nuclear  wave  function  which 


is  a  na.rrow  Gaussian.  We  obtain  then  one  tralectorv  (i.e,  the 
mean  trajectory)  and  this  is  in  agreement  with  our  classical 
prejudices  that  one  trajectory  should  describe  the  motion  of  one 
atom. 

However,  if  we  forget  that  we  have  ever  heard  of  the 
Newtonian  classical  mechanJ.cs  and  only  think  of  it  as  a 
computationally  advantageous  limit  of  quantum  mechanics  (which  is 
taken  by  using  narrow  Gatissian  wave  packets)  than  we  will 
generate  an  entirely  different  picture.  Since  quantonn  wave 
function  is  of  the  form 

il)(x^R;t)  »  x^{R)*j^(x;R)  + 

where  Xj^  and  Xj  st*"®  nuclear  wave  functions,  we  should  replace  Xj^ 
(R)  and  X2(R)  with  two  distinct  Gaussians.  By  forcing  them  to  be 
narrow  we  can  vise  time  dependent  Schrodinger  equation  to  generate 
two  classically  looking  equations  for  the  centers  of  the  two 
Gaussians.  We  have  now  two  trai jectories,  one  for  each  electronic 
state.  The  Guassian  G^  associated  with  the  electronic  state 
moves  on  a  potential  that  resembles  and  G^  on  one  which 

resembles  ^22'  However,  the  two  trajectories  are  coupled  by  a. 
potential  which  is  large  whenever  the  two  Gaussians  overlap  and 
whose  magnitude  depends  on  the  coupling  H^^. 

The  picture  emerging  from  such  a  description  is  more 

complex  and  more  successful  in  the  sense  that"  it  provides'  a 

classical  like  description  which  contains  most  of  the  important  • 

quantum  effects.  The  occupation  probability  is  given  by  the 

2 

amplitudes  of  the  Gaussians  and  if  P.  is  large  then  (G.)  is 

2  ^  ^ 
large  and  (G^)  is  small.  Asymptotically  the  two  packets  emerge 

with  two  different  momentum  distributions:  the  one  moving  on  the 

upper  state  has  a  smaller  mean  momentum  than  the  one  moving  on 

the  bottom  state.  This  reflects  the  fact  the  particles  emerging 

on  one  electronic  state  have  a  different  kinetic  energy  than 


those  emerging  on  the  other.  MTA,  on  the  other  hand  gives  a 
kinetic  energy  which  is  the  average  of  this  distribution. 

Furthermore,  the.  trajectory  running  on  state  1  and  that 
running  on  state  2  have  different  turning  points,  which  better 
mimic  the  behavior  of  the  two  nuclear  wave  functions.  MTA  has 
one  turning  point  which  is  an  average  of  that  of  the  two 
trajectories. 

A  more  striking  demonstration  of  the  need  for  double 
trajectories  can  be  made  when  we  include  a  decay  mechanism  by 
which  the  particle  can  lose  energy  by  exciting  phonons  or 
electrons  in  the  solid.  The  trajectory  running  on  can  be 

trapped  (forever)  it  It  loses  energy  to  excite  the  solid;  the  one 
running  on  cannot  be  trapped  since  the  state  2  is  repulsive. 

Thus  the  behavior  of  the  two  trajectories  is  extremely  different. 
If  a  mean  trajectory  is  used  for  such  situations  it  will  either 
get  trapped  or  it  will  escape  since  it  will  give  the  average 
behavior  of  the  two  trajectories.  If  the  amplitude  of  the 
escaping  Gaussian  is  the  larger  one  than  the  mean  trajectory  will 
escape;  in  the  opposite  case  it  will  get  trapped.  It  is  unlikely 
that  the  average  behavior  will  mimic  faithfully  the  more  correct, 
"bimodal"  behavior  of  the  two  trajectories  model.  How  well  MTA 
compares  to  a  double  trajectory  theory  can  be  decided  only  by  a 
detailed  numerical  investigation,  which  is  now  in  progress.^ 
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FIGURE  CAPTIONS 


The  diabatic  and  adiabatic  energy  surfaces  and  the  coupling 
used  in  the  present  calculations.  For  the  definitions  see 
Eqs.  (2.4)  and  (■2.7)-(2.9)  . 

The  mean  trajectories  (atom-surface  distance  as  a  function 
of  time)  for  various  values  of  the  Massey  parameter  X 
(defined  by  Eq.  .(2.5))  or,  equivalently,  Incident  kinetic 
energy  K  (in  eV) .  The  arrows  mark  the  turning  points.  The 
intersection  with  the  curve  R  -  5  gives  the  times  (marked 
on  the  figure)  when  the  particle  passes  through  the 
crossing  point  of  the  diabatic  curves.  The  turning  points 
are  the  points  where  R  changes  direction;  the  turning 
distances  are  given  near  the  horizontal  arrows  and  the 
turning  times  near  the  vertical  ones. 

The  effective  potentials  (full  lines)  for  various 

values  of  the  Massey  parameter  X  as  a  function  of  the 
surface-atom  distance.  The  arrows  mark  the  direction  of 
particle  motion.  The  adiabatic  energy  surfaces,  and 

^22  given  by  the  broken  lines. 
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The  probability  of  remaining  on  the  adiabatic  state  4^ 
as  a  function  of  time,  for  various  values  of  the  Massey 
paramter  X. 

The  asymptotic  value  of  the  probability  that  the  system 
starts  in  the  adiabatic  state  1  and  emerges  in  the 
adiabatic  state  2,  as  a  function  of  X.  The  solid  line 
gives  the  MTA  results  and  the  broken  line  those  of  the 
‘Nikitin  theory. 

The  same  as  Fig.  2  except  that  the  system  starts  in  the 
adiabatic  state  2 . 


The  same  as  Fig.  3  except  that  the  system  starts  in  state 


The  same  as  Fig.  4  except  that  the  system  starts  in  state 


The  probabilities  P®  that  the  system  starts  and  stays  on 
state  2,  as  functions  of  time.  The  full  line  gives  the 
results  of  the  trajectory  approximation  with  the  adiabatic 
potential  ^^2’  broken  line  give  the  mean  trajectory 

approximation  (MTA)  results. 

The  transition  probabilities  ®  function  of  the 

Massey  parameter  X  for  several  coupling  strengths  The 

full  lines  give  the  results  of  Kikitin’s  theory,  and  the 
broken  lines  give  the  MTA  results.  Xj  and  X^,  are  def.ined 
in  the  text. 


The  asymptotic  values  of  probability  P“  (the  probability 

«  A 

that  the  particle  starts  on  state  2  and  emerges  on  state  1) 
as  a  function  of  X  for  ^  -  0.3  given  by  the  mean  trajectory 
approximation  (full  line),  the  trajectory  approximation 
with  the  adiabatic  potential  (broken  line)  and  the 

Nikitin’s  theory  (chain  line). 


The  same  as  Fig.  11  except  p  «  0.15. 


The  maximum  asymptotic  values  of  probability  P^  as 
functions  of  X  for  various  values  of  p  (full  lines)  and  the 
value  of  (2SX^)~^  (broken  line). 


The  probability  P*  (that  the  particle  starts  and  remains  on 
state  2)  as  a  function  of  time  for  several  values  of  X 
(solid  lines),  as  computed  by  .MTA.  The  broken  line 
indicates  the  values  of  P®  below  which  the  condition  (2.24) 
is  satisfied. 


The  mean  trajectories  for  several  values  of  X  (full  lines). 
The  broken  line  indicates  the  crossing  point. 

The  effective  potential  for  X  -  30.02  (full  line)  and  the 
adiabafic  potentials  (broken  lines).  The  numbers  (1)  - 
(5)  indicates  the  time  sequence  in  which  the  particle 
moves  on  the  effective  potential  (see  the  text). 

The  asymptotic  values  of  the  probability  (that  the 
particle  starts  on  state  and  emerges  on  state  1)  as  a 
function  of  X  given  by  the  mean  trajectory  approximation 
(full  curve  and  circles)  and  by  Nikitin's  theory  (-•-•-). 
The  numbers  near  the  circles  indicate  how  many  -times  the 
trajectory  turned  back  towards  the  wall.  The  broken  lines 
show  the  results  for  at  the  time  when  the  particle  hits 
the  turning  point  that  will  send  it  back  towards  the 
surface  for  a  second  collision  (see  text). 

Same  as  Fig.  15a. 

Same  as  Fig.  15a. 

Same  as  Fig.  16a.  The  arrow  indicates  the  value  of  X  for 

2  z. 

which  1/2  -  de/2  -  ^22^^"^p^  that  the  particle  turns 

back  without  the  collision  with  the  wall. 


Table  1. .  The  values  of  the  parameters  used  in  the  presenr^  calculations 


U  “  0.025  Hartree  0.68  eV 

o 

At  =  0.005  Hartree  0.136  eV 

a  “0.4  Bohr  radius”^  0.756  A~* 

a~^  “  2.5  Bohr  radius  1.32  A 

R  “5.0  Bohr  radius  2.64  A 

o 

0  -  0.1,  0.15,  0.3  .  . 


m  “  42300  (in  the  atomic  unit) 

5  “  h^a*/mAe  “  7.56  x  lO”** 

__  ^ 

Table  2.  The  relationship  between  the  values /\and  the  incident  kinetic  energy 

X  Incident  Kinetic  Energy 


(Hartree)  (eV) 

2  0.825  22.4 

0 

5  0.132  3.59 

10  0.033  0.8976 

20  0.00825  0.224 

30  0.00367  0.0999 

36  0.00255  0.0694 


mV*  /2  -  AE/2(SX*  “  3.30/X*  Hartree  -  89.a/X^eV 
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